The locality of the square-root method for improved staggered quarks 
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We study the effects of improvement on the locahty of square- rooted staggered Dirac operators in 
lattice QCD simulations. We find the localisation lengths of the improved operators (fat7tad and 
ASQTAd) to be very similar to that of the one-link operator studied by Bunk et al., being at least 
the Compton wavelength of the lightest particle in the theory, even in the continuum limit. We 
conclude that improvement has no effect. We discuss the implications of this result for the locality 
of the n*'^-rooted fermion determinant used to reduce the number of sea quark flavours, and for 
possible staggered valence quark formulations. 

PACS numbers: ll.15.Ha, 12.38.Gc 



I. INTRODUCTION 

One of the biggest headaches in lattice QCD calcula- 
tions is the inclusion of quarks in a realistic and read- 
ily simulated manner. Staggered ferniions have the ad- 
vantage over Wilson-like quarks that their residual chi- 
ral symmetry protects the light quark mass from addi- 
tive renormalisation. Staggered fermions are also much 
cheaper to simulate than those from the overlap formu- 
lation. The disadvantage is that the formalism natu- 
rally yields Nf = 4 degenerate fermion species (known 
as "tastes"). In the continuum limit these tastes de- 
couple, but at finite lattice spacing, a, there are 0{a?) 
taste-changing interactions. Improvement of the action 
can systematically reduce these interactions, leading to 
smaller splittings in hadron mass multiplets 1|, a more 
physical Dirac spectrum J3, Q and the correct chiral sup- 
pression of the topological susceptibihty _4, ^j . 

True QCD, however, requires a light "2 -I- l"-flavour 
sea quark mass degeneracy. The usual way of tackling 
this is to replace the fermion determinant in the parti- 
tion function by its n'^-root, with n either 2 or 4, to give 
reduced-taste staggered fermions |a, Q ■ This approach is 
inspired by the continuum factorisation into independent 
tastes. It is not certain, however, whether this process 
yields a theory that is sufficiently local that we are confi- 
dent of the universality in the continuum limit, and will 
arrive at the correct theory, QCD, as a ^ 0. 

To answer this we seek to prove at least the existence 
of a local operator whose determinant matches that of 
the reduced-taste sea quarks. Explicitly formulating this 
operator would also provide a consistent description of 
the valence quarks — currently lacking — and allow an 
unambiguous perturbative analysis of the theory. Atten- 
tion so far has focused on the simplest candidate, the 



,th 



rooted Dirac operator, as (det p) 



l/n _ 



deiip 



l/r. 



The square-root of the simplest, one- link staggered op- 
erator was studied in |y|. Free fermions were shown to 
be non-local both at finite a and in the continuum limit, 
with a localisation length of order the inverse of the quark 
mass, m. The limits of the eigenvalue spectrum placed a 
similar upper bound on the interacting theory. Numeri- 
cal measurements showed the operator had a continuum 



limit localisation length at least as large as the correla- 
tion length of the lightest confined state, the Goldstone 
pion. 

In this Brief Report we investigate to what extent the 
non-locality (or otherwise) is due to the taste-changing 
interactions. Using a similar methodology to [Sj, we com- 
pare the locahty of the fatTtad and ASQTAD improved 
staggered operators at various lattice spacings. We show 
the size of the taste-changing interactions has almost no 
effect, and all the square-rooted staggered operators are 
non-local in the continuum limit on a scale at least that 
of the largest Compton wavelength. We discuss the im- 
plications of this for possible reduced-taste formulations 
for valence and sea quarks in the context of the expected 
properties of the infrared Dirac spectrum. 



II. METHOD 

An exponential locality length, Hoc, can be defined as 
follows 8J. Given a candidate Dirac operator lp{x, y), we 
apply it to a point source £,{y) = Syz to define a wave- 
function il>{x) = J2y Pi^^y)^iy)- The operator is expo- 
nentially local if the locality function f{r) is bounded 
by 



/(r) = max {|f/'(a;)|} < exp 

\\x-z\\=r \ rioc 



(1) 



We require rjoc ^ in the continuum limit to guarantee 
universality, as has been seen for overlap fermions j2||. 

We measure Eqn. ^ for the one-link, fatTtad and 
ASQTAD staggered quark formulations [Ifl Um ■ The lat- 
tices studied are described in Table ^ We use peri- 
odic boundary conditions and study three lattice spac- 
ings, choosing the quark mass such that the pion (light- 
est singlet pseudoscalar mesoii) mass is fixed in physical 
units. For compatibility with ISl, we choose m such that 
rom^ ~ 1.3 (or 520 MeV) [Illlllll. The lattice size 
is fixed around 1.5 fm, or four Compton wavelengths of 
the pion. 

We use the same m for calculations with the improved 
fermion formulations. Whilst the different renormalisa- 
tions will change the value of the pion mass, Ref. [13 



TABLE I: The Wilson gauge action ensembles studied. The 
/3 = 6.0 configurations are from 19]. Tadpole improvement 
factors uo come from the mean plaquette |2Q|. We set the 
scale using rg [2H|. 
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0.8680 


0.017 


6.0 
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221 
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0.8778 


0.010 


6.2 


24* 


94 


7.380 (26) 


1.63 


0.8851 


0.007 



shows this effect is only slight, and does not affect the 
conclusions of this study. 

If M is the staggered Dirac operator, then M'^M de- 
fined on the even sites is an Hermitian, positive-definite 
matrix representing Nj = 4 tastes |T5, 16]. We study the 

locality of ^ = VmUI using Eqn. Q. In this, | • | is 
the SU(3) colour norm, and || • || the periodic Li ("taxi 
driver") norm jg. 

The square-root function is approximated using a 
Tchebyshev polynomial jl^l of order Upoiy — 500. The 
analytic upper bound on the truncation error is 0(10^^) 
[l8| . Before discussing the results, we must first be sure 
that the truncation errors in the polynomial do not affect 
our estimates of the locality function. 



Measuring /(r) for various 



^ 500, the error 



^poly 

should decrease exponentially, and so 

A(r,7ipoiy) = |/(r,npoiy) - /(r, 500)1 

— ^g-Bnpoly ^i _ g--B(500-npoiy)' 

~ Ae-^"^'"y for < Up^iy < 500 (2) 
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FIG. 1: Relative accuracy of the localisation function using 
the polynomial approximation to the square root of the one- 
link staggered action. Measurements are on 16^32 for the free 
theory and /3 = 6.0. 
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FIG. 2: Scaling of the effective localisation length for 
reduced-fiavour Nf = 2 ASQTAD fermions. 



and A{r),B{r) > 0. We plot A at various fixed r in 
Fig. ^ for free fermions and for a representative gauge 
background at (3 = 6.0. Extrapolating the linear region 
to ripoiy = 500, we estimate the truncation error to be 
0(10"®). This is an order of magnitude smaller than 
the minimum of f{r). As improved fermions should lie 
somewhere between the free and the one-link interacting 
case, we use ripoiy = 500 in all our analysis. 

We may ask why the actual error is so much smaller 
than the theoretical bound. A heuristic argument is that 
the errors on the Tchebyshev approximation are heav- 
ily concentrated at the lower end of the spectrum. The 
eigenvalues of the Dirac operator are, by contrast, rel- 
atively sparse here. Assuming the localisation function 
is not completely dominated by the infrared modes, it 
is reasonable that truncation errors towards the middle 
of the operator spectrum (which are many of orders of 
magnitude smaller) are more representative of the error 
on f{r). 



III. RESULTS 



Using the locality function, we define an effective lo- 
calisation length rf^^{r) at distance r by 



'^loc 



(r-i+i - Ti) 



log(/(r,)//(r,+i)) 



(3) 



for each pair of subsequent distances r^ . 

We plot this function for the square-rooted ASQTAD 
operator in Fig. [21 scaling everything in units of rg . At 
small distances /(r) is almost a power law (with expo- 
nent somewhere between —3 and —4). Comparing our 
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FIG. 3: A comparison of the effective locality length for 
various reduced-flavour Nf = 2 staggered operators at /3 = 
6.0. 



results with larger volumes aL ~ 2.1 fm, only results for 
r/ro < 3 are free of finite volume effects. In this range 
there is no clear evidence for a crossover to exponential 
locality, so localisation length estimates must be treated 
as lower bounds only. Nonetheless, these bounds do not 
scale to zero with the lattice spacing, but rather remain 
constant. It therefore seems highly likely that the actual 
localisation length also remains finite in physical units in 
the continuum limit. 

We conclude, then, that the square-rooted ASQTAD 
Dirac operator is severely non-local in the continuum 
limit. In Fig. we compare at fixed lattice spacing 
the three operators. The effective locality lengths are 
extremely similar. This was also the case at the other 
lattice spacings. It would therefore appear that the im- 
provement programme for staggered quarks has no effect 
on the locality of the square-rooted Dirac operator. 



IV. DISCUSSION 

Understanding how to consistently reduce the number 
of staggered tastes for the sea and valence quarks is vital 
to realistic QCD simulations. A proposed solution for the 
sea quarks is to take the n^'^-root of the fermion determi- 
nant. The locality (and legitimacy) of the reduced-taste 
theory is not clear, but may be answered by attempt- 
ing to find a local reduced-taste Dirac operator with the 
same determinant. Using this same operator to describe 
the valence quarks (rather than the Nf = 4 staggered 
operator used at present) would also avoid the poten- 



tial mode mismatches seen for other mixed valence/sea 
fermion combinations |22l |. 

We have examined a candidate two-flavour Dirac op- 
erator, inspired by the n'^'-root trick used to reduce the 
number of sea quark tastes. We have seen that the lo- 
cality of this square-rooted operator is as bad for im- 
proved staggered formulations (fatTtad and ASQTAd) 
as it was for the one-link case ^] : the locality length for 
5.8 ^ /3 ^ 6.2 shows no signs that Hoc vanishes in the 
continuum limit. By contrast, rjoc appears to be at least 
the Compton wavelength of the lightest particle. 

In many areas improved staggered quarks show a defi- 
nite improvement over the one-link case [llyjIJiH- That 
there is almost no commensurate change in the locality 
suggests that the n'^-rooted operator is not a useful way 
of formulating the valence quarks. It is therefore unlikely 
to be a good choice for studying the locality of the deter- 
minant; a better choice would be an operator with more 
physical properties. 

Given a sufficiently smooth gauge background of topo- 
logical charge Q, a lattice Dirac operator should have 
an infrared spectrum characteristic of the number of 
tastes/flavours it represents. The Atiyah-Singer Index 
Theorem predicts that there should be (at least) Nf \Q\ 
near-zero modes J3, yl 123] . The remaining low- lying 
modes should lie in near-degenerate iVj-plets, with any 
splitting due to taste-changing interactions |2, y| . The 
multiplets should follow a universal distribution given, 
for instance, by random matrix theory [3, UM- Taking 
the n*'^-root will reduce neither the number of zero modes 
nor the multiplicity of other eigenvalues by factors of n. 
In addition, eigenvalue ratios (As )/(A/") can no longer 
agree with the universal distribution. 

By analogy with the reduction of 16 naive fermion 
species to 4 ultra-local staggered tastes, it seems likely 
that a physical reduced-taste staggered Dirac opera- 
tor will, if it exists, be obtained by projection rather 
than rooting. Unlike naive to staggered, however, taste- 
breaking interactions prevent this being exact at finite 
lattice spacing, but improved staggered spectra already 
show much closer taste symmetry Q, Q . It is not clear 
whether such a reduced-taste operator would be local or 
even renormalisable. It is, however, important to study 
this to understand the systematic uncertainty of using 
the n*'^-root method in precision lattice QCD simula- 
tions. 
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